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The initial stages of the quasiparticle decay in a Fermi liquid are governed by a time-scale distinct from 
the scattering rates as derived from the Fermi golden rule approach. We show that the initial decay is non- 
exponential and that it is determined by the zeroth spectral moment of the electron self-energy. We analyzed 
numerically a number of approximations for the self-energy by comparing with exact configuration interaction 
calculations for small finite system with fragmented states. A numerically simple approach for computing the 
spreading of the quasiparticle states for large systems is devised. 
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With the recent spectacular advances in light sources and 
ultrafast spectroscopic methods it has become possible to 
trace the quantum dynamics of electronic systems down to 
the atto-second time scale ([1-3] and references therein). One 
of the prime goals is the understanding of the nature of the 
formation and decay of electronic states. In this respect, nu- 
merical methods are feasible for few-electron systems only. 
For condensed matter one has to resort to different concepts 
such as the Landau's theory of Fermi liquids [4] that de- 
scribes (low-energy) long-lived excitations as quasiparticles 
(QP). Due to residual interactions QPs decay in time gener- 
ally. How this decay proceeds in time t is known under certain 
conditions only: A QP decays exponentially as ^ cxp(— 7t) 
with the rate constant 7(e) ~ e^ jep [5]. Here e is the quasi- 
particle energy and ep is the Fermi energy. This, however, 
holds true only for times t ^ I/7. The importance of this re- 
striction is illustrated by the following: Let us assume an ex- 
ponential decay at all t, thus the QP peak appears in frequency 
with a Lorentzian shape. This means that at the QP energy e 
the spectral function behaves as Ai^{uj\ e) = ^ (^_22, 2 ■ The 
standard deviation (T^(e) of the spectral density given by such 
a functional form diverges. 



short and long-time behavior, i.e. 



a^e) 



duj{uj — e) Ai^{uj; e) — > 00. 



(1) 



Explicit calculations for a three-dimensional (3D) homoge- 
neous electron gas (HEG) show that this divergence is spuri- 
ous [6] and that the zero, first and second spectral moments 
are indeed finite. The convergence of the integral Eq. (1) is 
governed by the high-frequency behavior of the spectral func- 
tion. Thus, the short-time limit of the single particle Green's 
function (from which Aj^ derives) is of special interest. In 
this respect quantum-kinetics indicates a quadratic decay in 
time [7]. Experimentally, atto-second resolution of electronic 
states in condensed matter has already called for a careful in- 
spection of this issue [8]. 

Here we present a spectral function that exhibits the correct 



dt 



A(t;e)*A="e-^*. 
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(3) 



These equations are exact and can be obtained nonpertuba- 
tively from very general physical considerations [9]. An at- 
tempt with a similar goal has been undertaken in Ref. [9], the 
resulting spectral function, however, violates the sum rules 
and has a shape with the spectral moments finite at any or- 
der, at variance with Ref. [6]. The spectral function given 
in this work fulfill all sum rules and comply with the ex- 
act short and long time-limits. The key ingredients are the 
imaginary part of the on-shell electron self-energy 7(e) ~ 
ImS((x) = e; e) and the decay constant cr^(e) as expressed 
in terms of the zeroth spectral moment of the self-energy 
S(") = - J^ dw |ImS(a;; e)|. The approach provides a 
recipe to compute the short-time limit of the electron corre- 
lation function on the basis of many-particle perturbation the- 
ory. In particular, we demonstrate how the decay constant 
(Eq. 2) can be computed diagrammatically to any desired or- 
der in the interaction. Conceptually, the problem should be 
addressed by the quantum kinetic theory. However in this for- 
malism, analytic calculations of the initial stage of the quasi- 
particle are not available and numerical approaches rely on 
further approximations [10-12]. 

For the decaying part of the spectral function we make an 
ansatz 



A{t; e) = cxp -7(e) 



t^ 



t + T{e) 



(4) 



which obeys the two limiting cases [Eqs. (2,3)] with cr^(e) ~ 
2j/t. The analytic form of the Fourier transform of this func- 
tion (A(oj; ^) ~ -^ Jq cos[(liJ — e)t]A{t; e)dt) is not known, 
however it is possible to prove that in the frequency domain 
it has exactly three finite spectral moments in accordance 



with [6]. The large lo expansion reads 



A{u-e) 



6 



TT (cj — e) 



4^2 



as uj ■ 



A{lo\ e) is normaHzed, i.e., 0th spectral moment is one [13]. 

According to Altshuler et al. [14] the initial stage of the 
quasiparticle decay always involves a formation of the two- 
particle-one-hole state {2plh). The rate of the process is given 
by \he. first collision time 1/ti and is determined by the corre- 
sponding Coulomb matrix elements or, in other words, by the 
available phase-space (the energy and the momentum must be 
conserved). The phase-space also determines in a crucial way 
the subsequent stages of the QP decay, which results in the 
creation of an increasing number of particles and holes, form- 
ing either localized or delocalized states in a Fock space. In 
the latter scenario the exponential decay is established after 
many generations of particles and holes have emerged. From 
these very general arguments it is obvious that the exponential 
decay requires a certain time to develop, which in our theory 
is determined by the parameter T(e). This time parameter cer- 
tainly exceeds the first collision time (T(e) ^ ti) obtained 
from the golden rule arguments applied to the bare Coulomb 
interaction (at the initial stages the screening is not efficient). 
This indicates that the time T(e) cannot be obtained from ei- 
ther the bare nor the screened interaction and is distinct from 
the relaxation time at the large-time limit (l/^{e)). 

To obtain T(e) let us recall the relations between the n}^ or- 
der spectral moments M^"^ of the single-particle Green func- 
tion and that of the self-energy 1?^' [6]: 



m(") = I, 
s(") 






- [M(1)]2. 



(5) 
(6) 



Soo is the frequency independent real part of the self- 
energy [15]. These matrix relations directly follow from the 
Dyson equation, and can be obtained in any basis ([6] used a 
plane-waves basis). £ is a diagonal matrix with the elements 
given by the zeroth-order state energies. For finite systems 
Hartree-Fock basis states are appropriate. Writing the ma- 
trix of the spectral functions in terms of the imaginary part of 
the single-particle Green function (A(cj) = -|ImG((x))|) and 
likewise for the spectral function of the self-energy (S(a;) ~ 
-|ImS((x))|), and using the superconvergence theorem [16] 
the matrices are cast as frequency integrals: 



/oo 
da;cj"AH, 
-OO 
/OO 
da;S(u;). 
-OO 



n = 0...2, 



(7) 
(8) 



In HF basis Sqo is rather small and is proportional to the dif- 
ference of the direct and the exchange Coulomb energy com- 
puted with the Hartree-Fock and exact density matrix, i. e. 
related to the deviation of the natural occupations from 1 or 
0. Thus, by virtue of Eqs. (5,7) we arrive at the conclusion 
that the Hartree-Fock energies in the first approximation are 



given by the center of mass of the spectral function. Likewise, 
by using Eqs. (6,7) we establish a formula for the matrix of 
standard deviations [cf. Eq. (1)]: 



(T^ = 



did S(uj). 



(9) 



Formally an exact representation of this positively defined ma- 
trix can be written in terms of the six-point response func- 
tion [17] (Fig. la). Instead, we will compute the standard de- 
viations for finite systems using its factorizations. 

To determine the set-in time of the exponential decay for 
3D HEG we apply Eq. (9) in a plane-wave basis to the states 
close to the Fermi surface. Calculations of Vogt et al. [6] and 
Farid [18] show that (T^(e) has two contributions: i) A local, 
momentum independent and ii) A non-local, momentum de- 
pendent. Their sum remains finite and positive at the Fermi 
momentum (^^(ej?) > 0). This tells us that the set-in time for 
the exponential decay for quasiparticles in the vicinity of the 
Fermi surface behaves as 



r{e) 



27(e) 



eFcr2(e^)' 



(10) 



The prefactor in front of e^ can be obtained analytically [19]. 
We also note that for HEG the spectral function consist 
of a quasiparticle peak with the oscillator strength less the 
unity (coherent part) surrounded by the satellites (incoherent 
part) [20]. Our approach goes beyond a description of the 
quasiparticle peak: fine detailes of both parts are smeared out 
preserving, however, the particle number and having a correct 
asymptotic behavior Calculation of higher-order satellites for 
comparison requires inclusion of diagrams of higher-orders in 
the screened interaction W and is extremely computationally 
demanding ([21, 22]). Therefore, we verify the performance 
of different approximations for the self-energy by comparing 
its zeroth moments computed by the configuration interaction 
(CI) approach for n finite electron system. 

As a prototypical system we consider the widely studied 
Naj^ cluster [23-25]. For the current purpose it is advanta- 
geous for several reasons: i) it contains a small number of 
electrons making it accessible to full CI [26, 27], ii) it can be 
seen as a generalization of 3D HEG to a finite number of par- 
ticles (jellium model [28]). The cr^ matrix can be computed 
exactly by exact diagonalization of the many-body Hamilto- 
nian. We use an algorithm by Olsen et al. [29] based on the 
graphical unitary group approach [30] for the generation of 
the restricted active space (RAS) and full CI Hamiltonians. 
The calculations are performed for each spin multiplicity sep- 
arately using spin-adapted basis functions [31]. In terms of 
the matrix elements of the creation and annihilation operators 
the spectral moments are expressed as 



_D'^ 



M 



(«) _ 



- E (4: 



'XPX''+ ^(e-)"Y«Y', (11) 



where the summation is performed over the Hilbert space 
of the ionized states (dimension D^~^) and electron attached 
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FIG. 1. Diagrams for the electron self-energy, a) exact expression in terms of the 2plh response function in which the entrance and exit 
channels cannot be separated by cutting one fermion line; b) {p-h)-p factorization leading to the GW approximation; c) (p-p)-h factorization. 



states (dimension 13^+^). The matrix elements of electron 
creation (aj^,) and annihilation (aa) operators 

XS = (pN+l|4|0N), y„' = (gN-l|a„|ON), (12) 

and the transition energies e+ = E^^^ — Eq, e^ = Eq — 
E^^^ are computed from the CI many-body states. 

Approximations for the self-energy operator can be ob- 
tained from the factorization of the 2plh six-points func- 
tion [15, 17] (Fig. la). If the particle-hole (p-h) Green's func- 
tion is treated exactly we obtain the so-called GW approxima- 
tion [32] (Fig. lb). Alternatively, this approximation can be 
obtained from the ^[G, W] variational energy functional [33] 
expanded in terms of the dressed electron propagator G and 
the screened Coulomb interaction W. A single diagram of the 
first order has to be considered. Finally, one obtains the same 
functional form by neglecting the three-point vertex function 
F in Hedin's equations [34]. It should be noted that we do 
not perform the self-consistent solution of Hedin's equations; 
instead, we compute exactly G and W from the exact one- 
particle and particle-hole propagators. They are given by 
the Lehmann representations in terms of many-body electron 
states [35, 36]. In accordance with the spectral representation 
of the self-energy [26] we obtain for the energy-uncertainty; 

^a!^''^ =^J2Y. («Pl">("b/3) + {p-h)-h terms, (13) 

where we introduced a notation for the convolu- 
tion of the Coulomb matrix elements {a(3\jS) = 

J d{rir2)(j)airi)(l)p{ri)(j)j{r2)(j)sir2)/\ri - r2| with matrix 
elements of the creation (or annihilation) operators [Eq. (12)] 
and with the density matrix elements Q"^ = (nNlaj^aalON); 






(14) 



Analogically, one obtains an expression for the self-energy us- 
ing the {p-p)-h factorization (Fig. 1 c); 

^N + 2 £)N-1 

m q 

where similar to Eq. (14) we define the convolution of 
{aji\^S) with the matrix elements of the two creation (annihi- 
lation) operators P'^^^ = (nN + 2\a\a\\m): 

{aq\m) =J2J2 >";(«/3|7'5)P;?. (16) 

/3 7^ 



Eqs. (13,15) can be thought of as the Fermi golden rule ex- 
pressions. Since the delta-function ensuring the energy con- 
servation is not present here the whole expression has a di- 
mension of the energy squared. Starting from Eqs. (13,15) 
we further derive a series of simpler approximations. When 
the HF Green's function is used in place of G we obtain 
the so-called G'^W approximation. If, furthermore, the non- 
interacting excited states are used to compute Q"^ we obtain 
the G'^W'^ approximation with the spectral moment; 






(17) 



where n^ is the occupation of the single-particle state 7. 

For the Na,^ cluster we computed the spectral moments 
(Eqs. 11,13,15) by exact diagonalization of the many-body 
Hamiltonian. The exact spectral function (solid line. Fig. 2a) 
of the lowest valence state {e^^ — 9.786 eV) is fragmented 
(two major peaks) and has multiple satellites. Despite this 
fact the model spectral function centered at the HF energy 
approaches the exact one in a large range of energies. The 
exact energy-uncertainty from the first two spectral moments 
of the spectral function (Eq. 11) is compared with the ex- 
pressions resulting from the approximation for the self-energy 
(Eqs. 17,13,15) (Fig. 3). Corresponding self-energy spectral 
functions and weights are in Fig. 2(b,c). Generally, GW self 
energy yields results superior to other approximations. De- 
spite the fact that {p-p)-h factorization performs bad for S(a;) 
(large energy gap in the case of finite systems) the spectral 
moments are close to those of G^W^ approximation. 

Summarizing, we presented a form of the quasiparticle line- 
shape that reflects the correct short and long time-limits of 
the single-particle Green function and, thus, can be used to 
parametrize evolving in time electronic structure (e. g., at- 
tosecond time-resolved photoemission). The spectral function 
also explicitely enters a description of a number of static pro- 
cesses (e. g. core state x-ray photoemission). Thus, the ex- 
perimentally observed form of the Fermi edge singularity will 
be affected by both the finite life-time of the core state I/7 as 
well as by the corresponding set-in time r [13]. 

In the case of 3D HEG the spectral function describes both 
the coherent and incoherent parts. We show that the set-in 
time r(e) vanishes as e^ in the vicinity of the Fermi surface. 
In the case of finite systems the CI method enables us to com- 
pute the energy uncertainties using a number of diagrammatic 
approximations. Our simulations indicate that accurate results 
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FIG. 2. (Color online) CI calculations for the Nag cluster: a) The exact (solid) and model (dashed) spectral functions of lowest valence state. 
The dots denote the weights (squares of the matrix elements, Eq. (12)) of the exact spectral function; b) The spectral strength of the self-energy 
in {p-h)-p (red) and {p-p)-h (yellow) factorizations; c) The spectral function (Tr[S(ci;)]) of the exact (solid), [p-h)-p (dashed) and {p-p)-h 
(dash-dotted) self-energies. 



are obtained even by neglecting the three-particle vertex, how- 
ever, further verifications are needed for extended systems. 
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FIG. 3. (Color online) Energy-uncertainty for the Na,]" states. Shaded 
area denotes deviation of the approximate theories from the exact. 
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General properties 

The spectral function of a fermionic many-body quantum 
system is defined in its most general form as the overlap of a 
particle (hole) states: 

2^S{t-t'-e) = {i:{t;e)\i,(t'-e)). 

The states emerge by creating a particle for e > ei? or a hole 
with e < ep and are not the eigenstates of the system in gen- 
eral and, thus, decay in time. It is convenient to represent the 
spectral function as a product of a noninteracting oscillatory 
part and a decaying part A{t\ e): 



S{t-e) 



2-K 



-A{t-e). 



For the latter we make an ansatz 



A{t; e) = exp -7(e) 



t + r(e) 



(1) 



The Fourier transform of the spectral function can be writ- 
ten as follows: 



A{u-e) 



cos[{uj-e)t]A{t-e)dt. (2) 



For A{t; e) given by Eq. (1) the analytic form of the Fourier 
transform is not known, but can easily be obtained numeri- 
cally [Fig. (1)]. Its odd spectral moments are zero because of 
symmetry consideration A{t; e) ~ A{—t; e): 

/•OO 

M(2fc-i)(g) ^ / doj{Lj- e)2*=-M(a;; e) = 0. 



Thus Af(^)(e) — M^^^{e) — e = 0. The even spectral moments 
can analytically be obtained from the derivatives of A{t; e) at 
t = 0: 

M(2fc)(g)= f duj{LU-ef''A{Lu;e) = (-l)''l[iaA^^''\t;e). 

This leads to M(°)(e) ^ M(°'(e) = 1 (normalization condi- 
tion), and M(2)(£) ^ M(2)(e) - [M^{€)Y = ^. Higher 



spectral moments diverge because of the discontinuity of the 
derivative of A{t\ e) at t = 0. This competes the prove that in 
the frequency domain it has exactly three finite spectral mo- 
ments in accordance with [1]. 

The asymptotic large lu expansion can readily be obtained 
by integrating Eq. (2) by parts : 



A{Lu;e) 



6 
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TT (w - e)4r2 



as uj 



Fermi edge singularities 

The proposed spectral function can be applied to obtain 
more accurate descriptions of the x-ray photoemission (XPS) 
line shapes from core states. The theory describing the 
many-electron response in XPS originated in works of An- 
derson [2, 3], Mahan [4], Nozieres and DeDominicis [5]. It 
was shown that a creation of an infinite number of electron- 
hole pairs accompanies the photoemission event. The result- 
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FIG. 1. up A{ijo\ 0) for the model spectral function defined by Eq. (1) 
for the rate constant 7 = 1.0 and different values of the set-in time. 
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FIG. 2. a) Line-shapes for the x-ray photoemission from core states with finite core state broadenings 7 and different set-in times r. r = 
(solid line) corresponds to the Doniach and Sunjic (DS) line-shape, b) The effect of finite set-in time (dotted) can results in a XPS profile 
experimentally indistinguishable from the DS line-shape with a substantially reduced broadening. 



ing electron spectrum for an infinitely-long-lived core hole is 
given by t^^ i1„ , where < a < ^ is called the Anderson 
singularity index and can be given in terms of the scattering 
phase shifts (5; as 



a 



2^(2/ + l) 



h 



For brevity we set the origin of the energy scale at the no-loss 
position of the hole-state. 

This approach is oversimplified and a proper treatment re- 
quires a finite lifetime of the core holes to be taken into ac- 
count. Convolution of the singular line-shape with a simple 
decay function cxp(— 7^) yields the well-known Doniach and 
Sunjic[6] profile: 



y°(^) 



2r(i 



(ixP- + 7^) 2 



1 /, X ^ 

-ira + [1 — a) arctan — 

2 V7 



where F is the gamma function. 

It is clear that even this modification is insufficient as it is 
based on the wrong assumption of the hole's spectral function. 
Therefore, we propose a modified line shape based on the pre- 
sented spectral function with correct asymptotic behavior. 



Yico) 



lim 

5-J-+0 



cxp 



It' 



r+ t 



dt 



{S - it)°' 



(3) 



where infinitesimally small positive S was introduced to select 
a branch of the multivalued l/(it)" function and to facilitate 
the numerical integration (Fig. 2). 

A comparison of XPS line-shapes for different set-in times 
(Fig. 2, b) reveals a similarity between the effect of finite set- 
in time and that of the reduced broadening. In other words, 
while 7 leads to the broadening of spectral lines the presence 
of non-zero set-in time t has an opposite effect as evidenced 



by M(2)(€ 



27(e) 

r(e) 



This observation might have a strong 



impact on the interpretation of the experimental XPS spectra. 
We also note a different approach leading to a modified 
spectral line-shape, the assumption of the frequency depen- 
dent scattering phase as was recently implemented for low- 
dimensional systems by Mkhitaryan and Raikh [7] . 
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